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1.0 SUMMARY 


A computer program entitled PAIN (an acronymn for Propeller 
Aircraft interior Noise) has been developed to permit calculation 
of the sound levels in a model of a cabin of a propeller-driven 
airplane. The model fuselage is a cylinder with a structurally 
Integral floor. The cabin sidewall and floor are stiffened by 
ring frames, stringers and floor beams of arbitrary configura- 
tions. The cabin interior is covered with a trim (i.e., insulation 
with a lining) to increase the sidewall sound Isolation and 
provide absorption in the cabin. 

The propeller noise of concern is actually a series of tones that 
occur at the blade passage frequency and at integral multiples of 
that frequency (i.e., at its harmonics). The program permits the 
calculation of the space-average interior sound levels for the 
first ten (10) harmonics of a propeller (of any design) rotating 
alongside the fuselage. 

The input data required by the program include the mechanical and 
acoustical properties of the fuselage structure and sidewall trim. 
Also, the precise propeller noise signature must be defined on a 
grid that lies in the fuselage skin. The propeller data have to 
be generated with a propeller noise prediction program for 
example, the NASA Langley ANOPP program [1]). 
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2.0 INTRODUCTION 


The present program mechanizes an analytical model presented 
in Reference [2]. Details of the underlying theory on which 
this program is based are available In that report. 

The two primary elements of the interior noise model are the 
fuselage and the propeller (Pigs. 1 and 2). The fuselage 
consists of a cylinder stiffened by ring frames and stringers, 
and a floor that is structurally an Integral part of the 
fuselage. The cabin space is the volume above the floor. The 
Interior surface of the cabin (sidewall) is finished out with 
a trim consisting of insulation covered with a lining. The 
propeller rotates about an axis parallel to the centerline of 
the fuselage. The model can be used to predict the space 
average sound levels in the cabin space for each of the 
various harmonics of the propeller (up to a maximum of ten 
(10) harmonics). 

The present model works with the pressure time histories (sig- 
natures) as defined over the fuselage at a number of closely 
spaced points on a grid that lies in the fuselage skin 
(Fig. 3). The pressure signatures are input as Fourier series 
that define the amplitudes and phases of each of the harmonics 
of the propeller (at each location on the grid). 

The fuselage structural modes are developed for the case of a 
stiffened cylinder with a floor partition. The structural 
properties of the cylinder and floor are required as input 
data. The structural modes are described by the computed 
resonance frequencies and mode shapes, and by associated 
structural loss factors. The loss factors of the modes of a 
bare fuselage must be input and must come from measurements 
(or be estimated). When trim is installed on the sidewall, 
the structural losses increase due to the trim's presence 
against the sidewall and this is computed for the particular 
trim installation. ^ 



T rim Insulation 
& Lining 
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FIGURE 1. PROPELLER AIRCRAFT INTERIOR NOISE MODEL 


Prooeiler Plane 



y 


FIGURE 2. PROPELLER AND FUSELAGE SURFACE POINT GEOMETRY 
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FIGURE 3. GRID USED FOR PROPELLER NOISE PREDICTIONS 






The acoustic modes of the cabin space are similarly defined by 
the resonance frequencies, mode shapes, and loss factors. The 
acoustic loss factors must be input for a bare fuselage but 
are calculated when a cabin trim is installed. 

The effects introduced by the trim at a given frequency are 
computed using the values of the wave Impedance and complex 
acoustic wavenumber in the insulation, and the trim lining 
surface weight and loss factor. 
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3.0 MODELS 


One can begin to construct an Input card package once the 
basic size of the fuselage is known and the location of the 
propeller defined. The overall geometry must be fixed so that 
the grid that will lie in the fuselage skin (Pig. 3) can be 
located. The grid is used for specification of the propeller 
noise excitation. The propeller-induced pressure signature at 
each location on the grid is required, and a propeller noise 
prediction program such as the NASA Langley ANOPP program [ 1 ] 
will have to be used to make these predictions. Obviously, 
the propeller diameter, rpm, number of blades, blade shapes, 
aircraft forward speed, etc., will be needed for these calcu- 
lations. The program PAIN uses free-fleld data as input and a 
coefficient is applied to the free-fleld amplitudes in a crude 
attempt to model the surface reflections caused by the fuse- 
lage. If a propeller noise prediction program that takes into 
account the reflection from a rigid cylinder located nearby is 
used, the surface reflection corrections applied in the 
present program can be removed and the new type of propeller 
data can be input. Just as is done at present with the 
free-fleld data. 

After the propeller data are generated, the remainder of the 
input requirements are associated with the definition of the 
fuselage structure and the trim on the cabin sidewall. 

The following is a step-by-step procedure used to parameterize 
the Interior noise model: 
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1) Referring to Pig. 2, define the propeller hub position along- 
side the fuselage by the coordinates (r ,<l> ). The dlmen- 

p P 

Sion rp is the radial distance from the center of the fuse- 
lage cylinder to the axis of rotation of the propeller. 

4> defines the angular position relative to the bottom of the 
fuselage. The location of the propeller relative to the 
front of the cylinder is given by Zp. The basic fuselage 
configuration is fixed by three dimensions — the cylinder 
length L, radius a, and the position of the floor as defined 
by the angle 6^ (see Pigs. 1 and 3). 


2) The grid used for the propeller noise predictions must next 
be located in the skin. The grid is defined over one-quarter 
of the periphery of the cylinder. Three coordinates of 
concern are the so-called ANOPP coordinates (x]^, X 2 , 

X 3 ) in Pig. 3 . Por purposes of defining the coordinates of 
the grid points (k, Z), the coordinate x^ , whose origin is 
at the propeller hub, will always penetrate the fuselage 
perpendicularly and that point of penetration will establish 
the line (Z = 1). The coordinate positions are then given by 
the following relations for k = 1, and = 1, 10; 


(k,£) (Xp, x^, x^) 

where x^ = (rp-a) + a{l-cos[ ( J.-1) tt/18] } 
x^= -a sln[ ( £-1 )tt/18] 

X ^ = (iTa/l8) [k -k] 


kp is the k index of the grid point (k,Z) = (kp,l), that 
is, where the coordinate xj penetrates the fuselage. Por 
example, for k = I 6 as in Pig. 3, setting kp = 8 
roughly centers the grid longitudinally about the propeller 
plane. Note that the value of (Tra/l8)*[kp - 1] cannot 
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exceed z„. Similarly, the chosen value of should be 

P TTlo. X 

such that (na/l8) x [k^^^-1] < L. This will assure that all 

rnd, A 

grid points are on the cylinder. 

3) The propeller noise at each location (k, l) will be period- 
ic with period T^ = BPP~^ where BPF is the blade passage 
frequency 


BPF 


NB 


where N is the rpm and B is the number of blades. The pro- 
peller noise at each location on the grid, over the period 

t = 0 to T]^, is (Eq. (33) of Ref. 2). 

^ (k , Jl ) °° 

+ n (k ,i!) ^ _ ^(k,£)^ 


p[ (k ,2.) ,t] = 




H=1 


(“ 






a(k,^')/2 is just the mean pressure and is of no concern. 
However, the two coefficients ’ and must be 

defined. Fourier analysis of the computed pressure time 
histories over the blade passage period T^ is necessary to 
provide the first ten (10) harmonics. The coefficients 
and for each position (k,Jl) are given by Eqs. (31) and 
(32) of Ref. 2. 


4) The structural properties of the fuselage cylinder and 
floor are required next. These data consist of the 
following : 

a) The elastic moduli! E and Poisson’s ratios v of the 
materials composing the cylindrical shell and the 
floor plate. 
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b) The equivalent thicknesses t of the shell and of 
the plate. These are given by the formula 



where is the actual thickness of the skin of the 
shell (or the plate thickness), Ag is the cross- 
sectional area of the stringers on the shell (or 
the longitudinally running floor beams), d is the 
spacing between stringers (or floor beams), Ar Is 
the cross-sectional area of the ring frames on the 
shell (or the transverse floor supports), and H Is 
the spacing between rings (or floor supports). Both 
the shell and plate must have a constant thickness, 
and the cross-section is uniform along the length of 
the cylinder. 

c) The bending rigidities of the stringers (or longitu- 
dinal beams) and the ring frames (or floor supports) 
divided by the spacings, l.e. 

Vstrlnger ’ “ 

a Ifloor 
/support 

The area moments of inertia of the stiffeners I are 
to be computed about the centerline of the skin (or 
plate) to which they are attached. If the stiffeners 
are not identical or not equally spaced, typical or 
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average values should be used for the bending rigid! 
ties/unit width. 

d) The average surface mass m(kg/m^) of the cylinder 
wall (or floor plate), l.e. 


m = pt, 

where P is the density of the material and t is the 
equivalent thickness as defined in (b) above. 


5) Having defined the structural properties of the fuselage 
cylinder, the modes of the fuselage are calculated using 
displacement series (see [2], Appendix D). The shell 
displacement in mode r is of the form 

n* 

'P^(z,6) = sln^^^^ ^ cosnS (or (-l)^"^^slnne ) 


for symmetric or antisymmetric modes, and for the 
floor. 


4^p(z,x) = 


sin 


Mttz 


n* 

E 

n=0 


.pr 

'Mn 


COS' 


niTx 


(or 


sin^^) 


where M = number of longitudinal half-wavelengths (Max 15) 
n = number of shell circumferential waves (Max 14) 
or number of half-wavelengths In floor width 
(Max 8). 

n* * number of terms In displacement series for shell 
or floor 


The symmetric modes and antisymmetric modes are calculated 
separately, for the same range of values of M and n. The 
cylinder is assumed to be supported at both ends by shear 
diaphragms. 
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6) Next come the geometrical properties of the cabin space 

which are given entirely by the cavity length L^, 
radius a, and the floor angle 9^. It is not necessary 
that the cavity length be equal to the structure length L. 
A variable d is defined that is the difference in the z 
coordinates between the start of the cavity and the start 
of the structure and is positive if L < . Usually L 

will be equal to or greater than and d will be zero or 
negative. 

7) There are certain data that pertain to the trim installa- 
tion, namely 

a) The wave Impedance of the insulation 

W = iW|e^'J’, 

and the complex acoustic wavenumber 

y = a - ik. 


where 


k = 2n /X , 

m ’ 

and X is the acoustic wavelength in the material. 

b) The trim lining mass per unit area m^ and the flexural 
loss factor of the lining. 

The wave impedance is given by the amplitude |W| -in 
mks rayls and the phase 4 > in degrees. The attenuation 
constant a is input in dB/m and the acoustic wave- 
length X^ in meters. The data are, of course, frequency- 
dependent. They are read in as a table at one-third 
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octave band center frequencies. In the case of propeller 
noise prediction, the blade passage frequency and its 
harmonics will fall at frequencies set by the propeller 
speed and number of blades. In such a calculation, the 
frequency bands have no purpose other than defining the 
frequencies for the data required for the trim. The 
insulation properties needed at the frequencies of the tones 
(l.e. the BPP and its harmonics) are Interpolated from the 
input at the band center frequencies. Data for a standard 
aircraft glass fiber insulation are given in Appendix A of 
Ref .[2]. 

8) When the acoustic and structural modal densities are large 
enough, the transmission can be expected to be dominated by 
modal response close to the excitation frequency. High 
frequency formulations have been developed in Ref. [2] as an 
approximation in this case. These approximations result in a 
considerable saving in computation time. 

If the change from low frequency to high frequency formula- 
tion is desired at and above a given frequency, this fre- 
quency must be specified in Program PAIN. Three options are 
available for this. 

a. If the input frequency is zero, the program estimates a 
change over frequency based on the volume of the cavity. 
The frequency is the 1/3 octave band center frequency in 
which the following empirical value lies. 


f 


4.1 c 


ol 


w 


Hz 
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where V = cavity volume (m^) 

Col = Interior speed of sound (m/sec). 

b. If the input frequency is greater than the maximum 

1/3 octave band center frequency or propeller harmonic 
frequency, the low frequency formulation is used 
throughout. 

c. If any other frequency is input, the program uses the 
high frequency formulation at and above the first 

1/3 octave band center frequency in which the input 
frequency lies. 

It should be noted that the high frequency formulation (expected 
value) is calculated at all possible frequencies and output as 
supplementary information, although not necessarily valid at low 
frequencies. For tone transmission, a bandwidth of 7% has been 
used to define structural modes close to the propeller harmonic. 
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4.0 OVERALL PROGRAM DESCRIPTIONS 


The master (main) program PAIN that computes the interior 
noise levels uses data that are generated by four auxiliary 
programs. One of these programs calculates the acoustic modal 
characteristics of the cabin, two other programs the struc- 
tural modal characteristics of the fuselage, and a fourth the 
propeller noise field. 

The acoustic modes are calculated with the program CYL2D. The 
only input required by CYL2D is a single card to define the 
floor angle 9^. The acoustic modes are determined for a two- 
dimensional case (i.e. q = 0) and a radius a = 1 meter. The 
resonance frequencies and mode shapes are computed and the 
mode shapes normalized for use by the main program. The 
values that the eigenvectors take on the periphery of the 
cabin are extracted and then ordered in sequence of boundary 
points. A code is used to identify symmetric and anti- 
symmetric modes (1 = symmetric, 0 * anti-symmetric). These 
conditioned data are written on a file (tape) for recall by 
the main program. Appendix C of Ref. [2] may be consulted for 
more details of the calculations performed by CYL2D. 

The structural modal data are provided by a program called 
MRPMOD whose output is also to a file (tape) read by PAIN. 

MRPMOD requires input data in the form of cards that define 
the number of generalized coordinates used in the structural 
analysis program, the maximum number of axial halfwaves 
(presently M = 10), the number of generalized coordinates 
retained to develop an accurate structural mode shape, the 
cylinder length, and finally the cylinder and floor masses per 
unit of area. In addition, MRPMOD reads two files (tapes) which 
are output from the fundamental structural program MRP (MRP 
stands for M. R. Peterson who developed the original analysis 
for the floor partitioned cylinder). The program MRP reads 
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data cards that define the properties of the shell and floor 
plate materials, masses per unit of area, stiffener bending 
rigidities, and alphameric Information that directs the 
program's calculations and output. MRP must be run 
twice — once for the symmetric modes and once for the anti- 
symmetric modes. The two files (tapes) created by these runs 
are read by MRPMOD. The files (tapes) contain the computed 
mode shapes and resonance frequencies. MRPMOD takes the 
output from MRP, calculates the generalized masses, rank 
orders the modes by frequency, and selects the most 
significant generalized coordinates, that Is, those needed to 
recreate the structural mode shapes where required In PAIN. 
Table D-1 of Appendix D of Ref. [2] gives an example of MRPMOD 
output. 

The propeller noise data must be calculated with a propeller 
noise prediction program, such as ANOPP [1]. This auxiliary 
program Is neither part nor parcel of PAIN as are the programs 
CYL2D, MRP, and MRPMOD. Propeller noise data must presently 
be Input to the main program as a card deck that contains (In 
a sequence to be specified exactly In Sec. 5) the location of 
the grid points (k,£) of Pig. 3 and the corresponding free- 
fleld Fourier amplitudes and phases of each propeller harmonic 
at each location of the grid. 

The programs CYL2D and MRPMOD do not call any subroutines. 

MRP and PAIN each has a large number of subroutines. The flow 
chart that follows In Figure 4 shows the basic sequence of 
Input manipulation, computation, and output In the main and 
auxiliary programs. 
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PROPELLER AIRCRAFT INTERIOR NOISE COMPUTER PROGRAMS 



FIGURE 4. FLOW CHART FOR PROGRAMS 














5.0 INPUT DATA 


The following Is a card-by-card description of the input for 
the programs CYL2D, MRP, MRPMOD, and the main program PAIN. 

All input is in metric units, with angles in degrees. 

5 . 1 Program CYL2D 

5.1.1 Input Data 

Card Columns Format Description 

1 1-10 P10.3 Floor Angle, 0 (see Figure 1) 

2 1 6/7/8/9 End of File 

The program CYL2D requires the IMSL library subroutines EIGRF, 
LINV2F, VMULFF, VSRTR and VSRTU. 

5.1.2 Output File, Tape 11 

The output file. Tape 11, must be cataloged or stored for input 
to program PAIN. Tape 11 contains 


Floor Angle 

Locations at which eigenvectors are specified 

First 20 eigenvalues (symmetric and antisymmetric combined) 

Generalized masses 

Eigenvectors for both interior and boundary locations 

5.1.3 Printed Output 

Eigenvector locations, for cylinder interior and boundary. 
Eigenvalues for first 20 modes (symmetric and antisymmetric 
combined) for cylinder with unit radius of 1 meter. Graphical 
presentation of eigenvectors. 

Generalized masses and eigenvectors. 
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5. 2 Program MRP 


5*2.1 Input Data 

Card Columns Format 

1 1-80 20A4 

2 1-33 A3 

3 1-10 AlO 

4 1-7 A7 

5 6-21 4A4 

36-43 2A4 

66-80 3A4,A3 

6 6-20 El 5.0 

7 6-10 15 

11-15 15 


Description 

Title cards (as many as required). 
These are printed directly on the 
output but not on Tape 7 output. 

END 

INPUT DATA 
GENERAL 

FREELY SUPPORTED 

Shear diaphragms at ends of cylinder. 
RIGID or HINGED 

Joint condition for Floor-Cylinder 
Junction. 

SYMMETRIC or ANTISYMMETRIC 
Symmetry about a vertical axial plane 
of the modes required. 

Angle 00 (degrees) at which the floor 
Joint Is located (e^ = 0 Is the bottom 
centerline). 

Starting value of M=1 In the dis- 
placement series. 

Final value of M (Maximum no. = 15) 
(Max.M X No. of eigenvectors < 225) 
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Card 

Columns 

Format 

Description 

7 

( Contd) 

16-20 

15 

Increment for M = 1. M = number of 
longitudinal half-wavelengths. 


21-25 

15 

Starting value of ng=0 In shell 
circumferential displacement series. 


26-30 

15 

Pinal value of ng for shell 
(Maximum no. = 14). 


31-35 

15 

Increment for ng = 1. 

n„ = no. of shell circumferential wave 


36-40 

15 

Starting value of np = 0 in floor 
width displacement series. 


41-45 

15 

Final value of np for floor (Maxi- 
mum no. = b) . 


46-50 

15 

Increment for np = l. 

n^ = no. of half-wavelengths in floor 
width. 

b 

6-10 

L5 

T = logical value true. Shell u and 
V degrees of freedom included. 


11-15 

L5 

Y = logical value true. Floor u and 

V degrees of freedom Included. 


l<^-20 

L5 

^ = logical value false. No extra 
printing required. 

9 

1-5 

A5 

SHELL. Shell structural data. 

10 

6-20 

E15.0 

Ex(N/m2) Modulus of elasticity in 
the axial direction. 
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Card Columns Format Description 


10 21-35 

(Contd) 


E15.0 E (N/m2) Modulus of elasticity In 

y 

the circumferential direction. 


36-50 

E15.0 

Poisson's ratio 

51-65 

E15.0 

Vg Poisson's ratio 

66-80 

E15.0 

^x6 Shear modulus. For 

Isotropic materials. If G^g Is blank; 
it is computed using G = E/2(l+v) 

6-20 

E15.0 

Thickness of shell (m). Includes 



"smeared-out" stiffener areas. 

21-35 

E15.0 

Shell mass/unit area (kg/m2) 
Including "smeared-out" stiffener 



masses . 

36-50 

E15.0 

Radius of shell (m) 

51-65 

El 5.0 

Length of shell (m) 

6-20 

E15.0 

D'^= ^ — (N.m) 

xs L ajgtringer I2(l-v^) 



Additional bending rigidity of shell 
stringers/stringer spacing. 

21-35 

E15-0 

r -1 E (t'-t 3) 

Dgr>= 5 ring" (N.m) 

frame I2(l-v^) 


Additional bending rigidity of shell 
frames/frame spacing. 

(Moments of inertia of stringers or 
frames to be computed about shell 
centerline) 


- 21 - 


Card Columns Format Description 


13 

1-5 

A5 

PLATE Floor structural data 

14 

6-20 

E15.0 

Ejj(N/m2) Modulus of elasticity in 




the axial direction 


21-35 

El 5.0 

Ey (N/m^) Modulus of elasticity in 




the transverse direction 


36-50 

E15.0 

Poisson's ratio 


51-65 

E15.0 

V Poisson's ratio 
•/ 


66-80 

E15.0 

®xy Shear modulus. If blank, 

it is computed using G = E/2(l+v). 

15 

6-20 

E15.0 

Thickness of floor (m) includes 
"smeared-out" stiffener areas. 


21-35 

E15.0 

Floor mass/unit area (kg/m2) in- 
cluding "smeared-out" stiffener masses 

16 

6-20 

E15. 0 

ml ~^i ^ 

Dvn= ^ long floor — (N.m) 

’‘P L'^Jbeam 12 (l-v>) 

Additional bending rigidity of longi- 




tudinal floor beams/floor beam spacing 


21-35 

E15.0 

d' = 1^1 floor ^ — (N.m) 

L support 12(l-v^) 


Additional bending rigidity of trans- 
verse floor support s/support spacing. 
(Moments of inertia of stiffeners are 
to be computed about the floor center- 
line) 
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Card 

Columns 

Format 

Description 

17 

1-3 

A3 

END 

End of Input Data 




The following cards control the oper 
atlons performed by the program. 

18 

1-17 

A10,A7 

GENERATE MATRICES 

19 

1-10 

AlO 

CONSTRAINT 

20 

1-5 

A5 

SHELL 

21 

1-5 

A5 

PLATE 

22 

1-3 

A3 

END 

23 

1-5 

A5 

SOLVE 

24 

1-5 

A5 

PUNCH Outputs to Tape 7* 

(Not Card Punch) 

25 

1-11 

All 

FREQUENCIES 

26 

1-12 

A12 

EIGENVECTORS 

27 

9-10 

12 

20 Number of eigenvectors 

(Generalized Coordinates) to be 
output (Maximum no. = 20, and 
Max.M X No. of eigenvectors < 225) 

28 

1-11 

All 

MODE SHAPES 
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Card Columns Format Description 


29 

1-4 

A4 

TRAN Transverse mode shapes 


38-40 

13 

20 Number of mode shapes to be 

output. 


43-45 

13 

3b Number of Increments for shell. 
[Displacements output at 5 degree 
Intervals] 


48-50 

13 

10 Number of Increments for floor 
[Displacements output at Intervals 
floor width Lp/20]. 


55 

LI 

T u and V displacements calculated 


60 

LI 

T output required on Tape 7. 

30 

1-3 

A3 

END 

31 

1-10 

AlO 

END OF JOB 

32 

1 

6/778/9 

End of file. 


No library subroutines are required. The program MRP must be run 
twice, once for SYMMETRIC modes and once for ANTISYMMETRIC modes. 
The same data Is used In both cases except for Card 5. 
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5.2,2 Output File, Tape 7 


The two output files. Tape 7, one for symmetplc and one for anti- 
symmetric modes, must both be cataloged or stored for input to 
program MRPMOD. These output files may be large if axial mode 
orders M = 1 to 10 are all included and care should be taken to 
allow for this. For example, on the GDC NOS computing service. 
Tape 7 should be defined as a Direct Access Pile rather than an 
Indirect Access Pile which is limited in size. 

For each value of M in turn, the output file. Tape 7, contains 
information on the first 20 normal modes including 

Resonance frequencies 
Generalized coordinates 
Shell displacements 
Floor displacements 

5.2.3 Printed Output 

Shell and Plate structural data. 

For each value of M in turn: 

Constraint matrices 

Eigenvalues and Resonance frequencies 
Generalized coordinates for first 20 modes 
Shell and floor displacements for first 20 modes 
Graphical representation of mode shapes. 
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5- 3*1 Input Data 


Card Columns Format Description 

1 l-80 16A5 Title card. This will be passed to 

program PAIN by Tape 9 output. 

2 1-5 15 No. of terms used for shell 

displacement series In program 
MRP [Final ng + 1] 

6-10 15 No. of terms used for floor 

displacement series In program 
MRP [Final np + 1] 

3 1-5 15 No. of values of M In program MRP. 

Maximum no. of longitudinal half- 
wavelengths. 

4 1-5 15 No. of terms used for shell dis- 

placement series In program PAIN 
[Maximum = 5] 

6-10 15 No. of terms used for floor displace- 

ment series In program PAIN 
[Maximum = 3] 

[For these, the most significant 
terms of the series will be selected 
by program MRPMOD] 

5 1-10 FIO.O Angle 0^ (degrees) at which the floor 

joint Is located (0^ = 0 at bottom 
centerline) 
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Card Columns Format Description 

6 1-10 PIO.O Cylinder length (m) 

11-20 FIO.O Cylinder radius (m) 

7 1-10 PIO.O Shell mass/unit area (kg/m^) includ- 

ing * smeared-out ' stiffener masses 
11-20 PIO.O Floor mass/unit area (kg/m^) includ- 
ing * smeared-out ' stiffener masses 

8 1-31 5X,I1,5A5 TAPE ? SYMMETRIC MODES 

9 1-31 5X,I1,5A5 TAPE 8 ANTISYMMETRIC MODES 

10 1 677/8/9 End of Pile. 

5.3.2 Input Plies 

Two files must be specified 

Tape 7 - Symmetric modes output from Program MRP 
Tape 8 - Antisymmetric modes output from Program MRP 

5 . 3.3 Output Pile, Tape 9 

The output file. Tape 9, must be cataloged or stored for Input to 
program PAIN. Tape 9 contains Information for all the modes, 
symmetric and antisymmetric for all values of M, arranged In 
ascending order of frequency. The output for each mode Is 

Frequency 
Mode symmetry 
Generalized mass 

Most significant generalized coordinates for shell and floor 
Shell displacements at 5 degree Intervals 
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5.3.4 Printed Output 


A summary is printed for each mode giving: 

Frequency 
Mode symmetry. 

M = number of longitudinal half-wavelengths 
C^n = Most significant shell generalized coordinates (Max 5) 

n = number of shell circumferential waves, associated with 
s 

each shell generalized coordinate 
C^n = Most significant floor generalized coordinates (Max 3) 
np = number of half-wavelengths in floor width associated with 
each floor generalized coordinate 
Total generalized mass for the mode (kg) 

Component of generalized mass due to shell normal displace- 
ment (kg) 

Component of generalized mass due to floor normal displace- 
ment (kg) 

5.4 Program PAIN 

5.4.1 Input Data 


Card 

Columns Format 

Description 

1 

1-80 16A5 

Title card 


2 

1-10 

F10.2 

Exterior 

air density (kg/m3) 


11-20 

F10.2 

Interior 

air density (kg/m3) 


21-30 

F10.2 

Exterior 

speed of sound (m/sec) 


31-40 

F10.2 

Interior 

speed of sound (m/sec) 

3 

1-5 

15 

NI = Number of 1/3 octave band 


frequencies (maximum 26) 
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cards will vary depending on the 
put format for frequencies, loss 
The card numbers that follow are 



Lowest 1 /3 octave band center 
frequency (Hz) 

2nd 1/3 octave band frequency (Hz) 
[8 frequencies per card] 


(Nl)th 1/3 octave band freq. (Hz) 

L = Fuselage cylinder length (m) 
a = cylinder radius (m) 

Bq = angle at which floor joint is 
located (degrees) = 0 is the 
bottom centerline 

V = Poisson's ratio (for end caps) 

[At the moment, transmission through 
the end caps has not been included 
in the program, but this may be added 
later] . 

Title card, describing the cavity 

Lq = Cavity length (m) 
d = distance of start of structure 
from start of cavity (d + ve L<L(^; 
d - ve L>Lq) 
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Card 

10 


1 1 

12,13, 


15 


16,17,1 


Columns Format Description 

1-10 F10.2 The center frequency of the first 

1/3 octave band in which the high- 
frequency formulation is used (Hz) 

( See Section 3 (8) ) 

If this value is zero , the program 
estimates a frequency. 

If this value is non-zero , the input 
value is used, not the estimated 
value . 


1-25 5A5 STRUCTURAL LOSS FACTORS 


4 1-10 F10.2 

11-20 P 10.2 


Structural loss factor for lowest 
frequency band. 

Structural loss factor for 2nd band. 


FI 0.2 Structural loss factor for (Nl)th 
band . 


1-25 5A5 ACOUSTIC LOSS FACTORS 

The acoustic loss factors must be 
input for a bare fuselage. If the 
input loss factors are zero, then the 
cabin trim data is used to calculate 
them . 


8 1-10 F10.2 


11-20 FI 0.2 


Acoustic loss factor for lowest 
frequency band (zero if calculated 
using trim data). 

Acoustic loss factor for 2nd band. 


F10.2 Acoustic loss factor for (Nl)th band. 
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Card 

Columns 

Format 

Description 

19 

1-80 

1 6A5 

Title card, describing the trim. 

20 

1-10 

P10.2 

Rt = insulation thickness (m) on 
curved wall of cylinder. 


1 1-20 

P10.2 

ni-t = mass/unit area (kg/m2) of 
the trim lining, facing the cavity 


21-30 

P10.2 

A-t = Surface area (m2) of the 
trim on the curved cylinder wall. 


o 

1 

P10.2 

= trim lining loss factor (in 


flexure) . 

***If there is no trim, the value of h^ = 0 must be input and 
Cards 21-33» the remaining trim data, should be omitted. 

With no trim, the acoustic loss factors must be input on cards 
16-18.*** 


If ht / 0 


21 


1-10 
1 1-20 


21-30 

31-40 


P10.2 he = insulation thickness on end 
caps (m) 

F10.2 mg = mass/unit area (kg/m^) of 

the end cap trim lining, facing the 
cavity . 

PI 0.2 Aq = surface area (m2) of the 
trim on one end cap. 

P10.2 rig = end cap trim lining loss factor 
(in flexure) 


**»Por each frequency band, the trim Insulation characteristics 
are input. (Por Piberglas Type PP105, the values are given in 
Ref. [3], Page 482)*** 
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Card Columns Format 

22,23,24 1-10 F10.2 

11-20 F10.2 

P1 0.2 

25,26,27 1-10 P10.2 

11-20 F10.2 

P10.2 

28,29,30 1-10 F10.2 

11-20 P 10.2 

PI 0.2 



a = attenuation constant (dB/m) for 
first frequency band. 


a for 2nd frequency band. 


a for (Nx)the frequency band. 

Xjjj = acoustic wavelength in material 
(m) for first frequency band. 

for 2nd frequency band. 


Xj^ for (Nl)th frequency band. 

I wl = amplitude of wave impedance 
(mks rayls) of trim insulation for 
first frequency band 

I Wj for 2nd frequency band. 


1 W| for (Nl)th frequency band. 
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Card 

Columns 

Format 

Description 


31,32,33 

1-10 

F10.2 

(j) = phase of wave Impedance (degrees) 
for first frequency band [+ve values] 

- 


11-20 

F10.2 

• 

j 

^ for 2nd frequency band. 




• 

F10.2 

(t> for (Nl)th frequency band. 


34 

1-15 

3A5 

PROPELLER TONES 


3^ 

1-80 

16a5 

Title card, describing the propeller 
data. Propeller location is 
(r*p,(}) ,zp) [See Figure 2] 


36 

1-10 

F10.2 

Tp = distance of propeller center- 
line from fuselage centerline (m) 



11-20 

F10.2 

(t> = angle between vertical and the 

line from propeller to fuselage 
centerline (degrees) [<{) must be given 
to the nearest 5 degrees] 



21-30 

FIO. 2 

Zp = distance of propeller plane 
from forward end of fuselage (m) 



31-40 

F10.2 

B = No. of blades of propeller. 



41-50 

FIO. 2 

Propeller rpm, for which propeller 





sound pressures were predicted. 
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Card 


Columns Format 


Description 


36 51-60 F10.2 +1.0 If propeller rotates counter- 

clockwise In z direction (as In 
Figure 3); or -1.0 if propeller 
rotates clockwise In z direction. 

(Input data created using ANOPP must 
be blade downsweep only) 


***The propeller pressure signature Is defined at grid points on 
the upper quarter of the fuselage skin, and Input, using a 
Fourier series analysis, as the amplitude and phase at the 
propeller harmonics for each grid point. The output from a 
propeller noise prediction program [1] has been used as part of 
this data Input package. The location of the grid points Is 
defined In Figure 3.*** 


37 1-5 15 


= No. of propeller harmonics 
max ^ ^ 

( Maximum = 10). 


6-10 15 k = Number of axial grid coordi- 

max 

nates (Maximum = 16). 


11-15 15 


kp = Axial grid coordinate at 
which propeller plane is located . 


The following x xIO) cards contain the pressure 
amplitude and phase for every grid point (k,^-) for all the 
propeller harmonics H. 


38 1-4 14 

onwards 5-8 I 4 


k = Grid axial coordinate 

2. = Grid circumferential coordinate 
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Card 

Columns 

Format 

Description 

38 

9-14 

F6.2 

xi' 

ANOPP coordinates 

onwards 

15-20 

P6.2 

X2 

, defined in Figure 3 

( contd) 

21-26 

P6.2 

X3 

(not used by Program PAIN) 


27-32 

16 

H = 

harmonic number. 


55-48 E16-5 <p^(k,X-) = phase of Fourier component 

for harmonic H at location (k ,il ) (degrees) 
( k Z ) 

49-64 EI 6.5 ~ amplitude of Fourier compon- 

ent for Harmonic H at location (k,ji) (Pa) 

This format is repeated until all values of k, Z and H have been 
used . 

1 ^/llQj^ End of file. 

No library subroutines are required. 


5*4*2 Input Piles 

Two files must be specified. 

Tape 9 - Structural modes output from program MRPMOD. 

Tape 11 - Acoustic modes output from program CYL2D. 

5 *4. 3 Printed Output 

1 . Band-averaged loss factors 

Acoustic loss factor , either input or calculated values when 
trim is present 

S 1 C 

structural loss factors input 

External radiation loss factors 

Total structural loss factor n^ when trim is present. 
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2 . 


Trim Data 


Trim characteristics input 

Trim transmission coefficient (dB) 10 log (1/t^) 

Trim and end cap admittances, real and imaginary. 

(non-d imensional) 

3. Number of structural and acoustic modes in each 1/3 octave 
frequency band. 

4. Structural Modes. For each mode: 

Frequency 

Symmetry 

Total structural loss factor 
M = number of longitudinal half-wavelengths 
For the most significant displacement series terms for 
shell and floor:- 

ng = number of circumferential wavelengths for shell 
^Mn = generalized coordinate for shell 
np = number of transverse half-wavelengths for floor 
^Mn = generalized coordinate for floor 
Total generalized mass (kg) 

Gen. mass associated with shell normal displacement (kg) 
Gen. mass associated with floor normal displacement(kg) 
Reverberant field joint acceptance j2(i’ev) 

5. Acoustic modes. For each mode: 

Frequency 

Symmetry 

q = no. of longitudinal half-wavelengths 

i = mode counter for cylinder cross-section modes. 

Acoustic loss factor n • 

n 

Normalization parameter 
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Acoustic/structural coupling factors. 


For acoustic mode n = (q,l) and structural mode r = (M,N) let 
f'(n,r) = fqM* fiM (see Volume I, Section 3.3) 

Values of f^^^ output for both shell and floor. 

Propeller Pressure Amplitude and Phase 

Output for all grid points, both above and below 6 = (|), for 

each propeller harmonic, 

Pressure amplitude (dB re 20 microPa) including reflecting 

surface effects. 

Pressure Phase (degrees) 

Propeller generalized forces. 

For each harmonic H and structural mode r. 

Generalized force 'i'„(r,H) - (see Ref. [2], Equation M2). 

U 

Tone Transmission for Propeller Harmonics - 

Interior space-averaged mean square pressure (dB re 20 microPa) 

1) Low frequency formulation calculations only, for all harmon- 
ics, identifying 5 modes making the largest contributions. 
Hp' output to identify closely coupled structural and 
acoustic modes. 

2) Summary of low frequency formulation (where specified) and 
the high frequency formulation (expected value) for all 
harmonics . 

3) Summary of predicted interior space-averaged mean square 
pressure (dB re 20 microPa). 


6.0 CONTROL CARDS FOR EXECUTING THE PROGRAM 

The control cards that follow are for the CDC CYBERNET NOS 
Computing Service. All the programs, CYL2D, MRP, MRPMOD and 
PAIN, are written in Fortran Extended 4, and comply with ANSI 
1966 Fortran for the most part, except where specific CDC 
functions such as EOF were required. In addition, it is 
assumed that IMSL library subroutines are available for 
program CYL2DF. 

The control cards assume that the Fortran programs are stored 
as permanent files on disc storage with the names CYL2DF, 

MRPF, MRPMODF and PAIN. Input and output data are also stored 
as permanent files. Sample input data for each program are 
also given. 
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ORlGtNAL PAG^ 

QE POOR QUALITY 


6.1 Control Cards and Data for Program CYL2D 


/JOB 

PAIN. T150.P4. 

/USER 

/CHARGE 

ROUTE, OUTPUT, DC=PR. UN=»******, ST=»**, DEF. 

GEI,CYL2DF . Fortran Progj’ain CYL2DF 

FTN,R=3. 

REWIND. 1 GO. 

RFPLACE,LGn=CY1.2DC. Com{-Ued CYL2DC 

GOTO, 1 . 

EXIT. 

1 , SI IME . 

DA YE I EE . 

/EOF 

/READ, CYL2DF 
/EOF 


/JOB 

PAIN. 1400. P4. 

/USER 

/CHARGE 

ROUTE . OUTPUT , DC=PR, UN=******* . ST=**« . DEF . 


GET,CYL2DD. 

ATTACH. IMSI /UN=LIBRARY. 

IBRARY, IMSL . 

GET. CYL2DC. 

COPYBF , CYL2DC, EGO. 

I DSEl ( PRESET = ZERO) 

I GO. PI =SOOOO. 

REPL ACE, TAPEl 1=CYL56. 
GOTO, 1 . 

EXIT. 

1 , STIME. 

DAYFIL E. 

/FOR 

/READ, CYl 2DD 
/EOF 


Input Data File CYL2DD 
Library Subroutines 

(-ornpl](?d Lronriui! CYIGIX: 
Output file CYLS6 


CYL2DD 


Input Data Pile 


56 . 6 
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6.2 Control Cards and Data for Prn p ;r»ain mrp 


/JOB 


PAIN.T150,PA. 

/USER 

/CHARGE 

ROUTE, OUTPUT, DC=PR, UN==******* 
GET.MRPF. 

FTN,R=3. 

REWIND, LGO. 

REPL ACE, LGO=MRPC. 

GOTO, 1. 

EXIT. 

1 , STIMF. 

DAYFILE. 

/EOR 

/READ,MRPF 

/EOF 


ST=***, DEF. 

Fortran Program MRPF 

Compiled Program MRPC 


Control Cards for Symmetric Modes 

/JOB 

PAIN, T400,P1 . 

/USER 

/CHARGE 

ROUTE. OUTPUT, DC -PR, UN=ihhhmhm(, ST = ***, DEF. 

DEF INE . TAPE7=MRPSM. Direct Access Output File MRPSM 

GET,MRPSD. Tnriut. Data File MRPSD 

GET, MRPC. Cornfjiiled Prograiri MPH,' 

COPYBF,MRPr.,l GO. 

L DSET ( PRESF T = ZERO ) 

LGO, PI =50000. 

GOTO, 1 . 

EXIT. 

l.STIME. 

DAYFILE. 

/EOR 

/READ, MRPSD 
/EOF 
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Control Cards for Antisymmetric Nodes 


ORJGSf^AL F/U-: !■.; 

OF POOR QUALITY 


/JOB 

PAIN,T400,P1 . 

/USER 

/CHARGE 


ROUTE, OUTPUT , DC=PR, UN=***#***, 
DEFINE. TAPE7=MRPAM. 

GET.MRPAD. 

GET.MRPC. 

COPY BE , MRPC , LGO. 

LDSET (PRESET=ZERO> 


S7=***. DEF. 

Direct Accecs Output Pile MRPAM 
Input Data File MRPAD 

Compiled Program IVRPC 


LGO, PI =50000. 
GOTO, 1 . 

EXIT. 

1 , STIME. 
DAYFILE. 

/EOR 

/READ, MRPAD 
/EOF 
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MRPSD Input Data File for Symmetric Modes 


STIFFENED .032 IN CYLINDER. 20 IN RADIUS, WITH FLOOR AT 56.6 DEGREES. 

MODE! #5 CYLINDER. PHASE III. 

FREELY SUPPORTED END CONDITIONS, RIGID JOINI AT FLOOR-CYLINDER JUNCIION 
SYMMETRIC MODES. OUTPUT ON MRPSM FOR M=^l-10 MODES 


END 

INPUT DATA 
GENERAL 

FREELY SUPPORTED 
56.6 

1 1 0 1 

T T F 

SHELL 

7. 239^95E10 
0.001 153058 
62. 352938 

PLATE 

7. 239'195E10 
O. 00127622 
8.6297699E3 


END 

GENERATE MATRICES 

CONSTRAINT 

SHELL 

PLATE 

END 

SOLVE 

PUNCH 

FREQUENCIES 
EIGE NVECTORS 


RIGID 

O 12 1 O 


7.239495E10 
3. 1 13933 
1 . 2O0725E4 


7.239495E10 
3.4465372 
5. 4264354E3 


1 


0 . 508 


SYMMETRIC 


1 . 8034 


15 


MODE SHAPES 

TRAN 15 36 10 T T 

END 

END OF JOB 
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MRPAD Input Data File for Antisymmetric Modes 


STIFFENED .032 IN CYLINDER, 20 IN RADIUS, WITH FLOOR AT 56.6 DEBREES. 

MODEL #5 CYLINDER. PHASE III. 

FREELY SUPPORTED END CONDITIONS, RIGID JOINT AT FLOOR-CYL I NDER JUNCTION 
ANTISYMMETRIC MODES. OUTPUT ON MRPAM FOR M=l-10 MODES 

END 

INPUT DATA 
GENE-RAL 

FREELY SUPPORTED RIGID ANT 1 S YMMt T R I ( 

56 . 6 

1 1<> 1 O 12 1 O 5 1 

T T F 

SHELL 



7. 239495E10 

7. 239495E10 


0. OOl 1 53058 

3.11 3933 


62. 352938 

1 . 2007 25E 4 

PLATE 




7, 239495F10 

7. 239495E10 


0. 00127622 

3.4465372 


8. 6297699E3 

5. 4264354E3 

END 



GENERATE 

MATRICES 



CONSTRAINT 

SHEL I 

PLATE 

END 

SOI VE 

PUNCH 

FREQUENCIES 

EIGENVECTORS 


0 . 508 


1 . 8034 


1 5 

MODE SHAPES 

TRAN 15 36 lO T T 

END 

END OF JOB 
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6 . 3 Control Cards and Data for Program MRPMQD 


/JOB 

PAIN, T 150, Pfl. 

/USER 

/CHARGE 

ROUTE , OUTPUT , DC=PR , UN=******» , ST=*** , DEF . 
GET,MRPMODF. Fortran Program MRPT'40DF 

FTN, R=3. 

REWIND, LGO. 

REPL ACE,LGO=MRPMODC. Cornp.Ued Program MRPMODC 

GOTO, 1 . 

EXIT. 

1 ,STIME. 

DAYFILE. 

/EOR 

/READ, MRPMODF 
/EOF 


/JOB 

PAIN, T400,P3. 

/USER 

/CHARGE 


ROUTE, OUTPUT, DC=PR, UN=*******, ST=***, DEF. 


GET.MRPMODD. 

ATTACH. TAPE7=MRPSM. 
ATTACH, TAPEa=MRPAM. 
GET.MRpMonr:. 

COPYBF , MRPMODC , LGO. 
L DSET (PRESET = ZERO) 


Input data file MRPMODD 
Cymmetrlc Modes File MRPFM 
Ant j syrnmctr‘1 c Modes P'lle MRPAM 
(;cjrn[) i ] cd Prajg.r'am MRPMODC 


LGO, PL =50000. 

REPL ACE.TAPE9=STR56. Output File STR56 

GOTO, 1 . 

EXIT. 

1 , STIME. 

DAYFII E. 

/EOR 

/READ. MRPMODD 
/EOF 


MRPMODD Input Data Pile 


STIFFENED .032 IN CYLINDER, WITH STIFFENED FLOOR AT 56.6 DEGREES 
13 6 

lO 


56.6 

1.8034 0.508 

3.11324 3.44578 

TAPE 7 SYMMETRIC MODES 
TAPE a ANTISYMMETRIC MODES 


ORlG^riAL PAGE 53 

OF POOR QUALITY 
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ORIGINAL FAOS 13 

OF POOR QUALITY 


6.^ Control Cards and Data for Program PAIN 


/JOB 

PAIN, T150, P4. 

/USER 

/CHARBE 

ROUTE, OUTPUT, DC=PR,UN=***#»**,ST=***. DEE. 

GET, PAIN. Fori. ran Program I’MIJ 

FTN, R=3. 

REWIND, L GO. 

REPLACE, LGO=PA INC. Compiled Program PAINC 

GOTO, 1 . 

EXIT. 

1 ,STIME. 

DAYFILE. 

/EOR 

/RTAD, PAIN 
/EOF 


/JOB 

PA I N , T 400 , P3 , CM370000 . 

/USER 

/CHARGE 

ROUTE, OUTPUT, DC=PR, UN=***»*** . ST=*** . DEF. 


GET, PAIND. 

GET, TAPE9=STR56. 
GET, TAPE 1 1=CYL56. 
GET, PAINC. 

COPYBF , PA I NC , LGO. 
LDSET (PRESET=ZERO) 
L GO , PL =50000 . 

GOTO, 1 . 

EXIT. 

1 . STIME. 

DAYFILE. 

/EOR 

/READ, PAIND 
/EOF 


Jrifjut Data FiJo I’AIHI) 
structural Modes File STR56 
Acoustic Modes File CYL56 
Compiled Program PAINC 
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6.M Control Cards and Data for Program PAIN 


/JOB 

PAIN,T150,P4. 

/USER 

/CHARGE 

ROUTE, OUTPUT, DC=PR, UN=*******, ST=***, DEF. 

GET, PAIN. Fortran Propxajn PAIN 

FTN,R=3. 

REWIND, LGO, 

REPLACE, LGO=PA INC. Compiled Program PAINC 

GOTO, 1 . 

EXIT. 

l.STIME. 

DAYFILE. 

/EOR 

/READ, PAIN 
/EOF 


/JOB 

PAIN, T400, P3, CM370000. 

/USER 

/CHARGE 


ROUTE, OUTPUT, DC=PR 
GET,PAIND. 

GET, TAPE9=STR56. 
GET,TAPE11=CYL56. 
GET, PAINC. 

COP YBF , PA I NC , LGO . 
LDSET (PRESET=ZERO) 


, UN=******#, ST=***, DEF. 

Input Data Pile PAIND 
Structural Modes Pile STR56 
Acoustic Modes File CYL56 
Compiled Program PAINC 


LGO, PL =50000. 
GOTO, 1. 

EXIT. 

1 ,STIME. 
DAYFILE. 

/EOR 

/READ, PAIND 
/EOF 
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APPENDIX A 


LIST OP SYMBOLS 

Ag (m^) Surface area of trim on one end cap. 

( k a) 

’ (Pa) Amplitude of Fourier component of propeller 

pressure signature at propeller harmonic H and 
grid location (k,^-) 

At, (m^ ) Cross-sectional area of ring frames on the shell 

A 

(or transverse floor supports) 

Ag (m^ ) Cross-sectional area of stringers on the shell (or 

longitudinal floor beams) 

A, (m^) Surface area of trim on the curved cylinder wall. 

L/ 

a (m) Cylinder radius 

k Z 

Sq* /2(Pa) Mean propeller pressure amplitude at grid location 
(k,£) 

B Number of propeller blades 

BPF = NB/60 (Hz) Blade passage frequency 

Cm^.c£^ Floor and shell generalized coordinates for 

Mn’ Mn 

structure mode r - (M,N) 

Cqi (m/sec) Interior speed of sound 

d (m) Spacing between stringers (or floor beams) 

d (m) Distance of start of structure from start of 

cavity [d + ve L < L^,; d - ve L > L^] 


A-1 


LIST OP SYMBOLS 
(Continued) 


E (N/m^) 
(N/m') 

Ey (N/m') 

Eq (N/m') 

f ' (n, r) = 

(N/m^ ) 
H 

IJ 

max 

hg (m) 

I (m‘) 

1 

j ^ (rev) 


Stiffener modulus of elasticity 

Shell or floor modulus of elasticity in axial 
direction 

Floor modulus of elasticity in the transverse 
direction 

Shell modulus of elasticity in circumferential 
direction 

Acoustic/structural coupling factor for 
acoustic mode n = (q,i) and structural mode r = 
(M.N) [see Ref. [2] Section 3-5] 

Shell or floor shear modulus 

Propeller harmonic order 

Maximum number of propeller harmonics considered 

Insulation thickness on end caps 

Insulation thickness on curved wall of cylinder 

Moment of Inertia of stiffener about shell or 
floor middle surface 

Acoustic mode number counter for cylinder cross- 
section modes, associated with mode n = (q,l) 

Reverberant field joint acceptance for each 
structural mode 


A-2 



LIST OP SYMBOLS 
(Continued) 


Axial grid coordinate (non-dimensional) for point 




Axial grid coordinate at which propeller plane is 
located 


max 


Maximum number of axial grid coordinates 


k = 2tt/A Acoustic wave number 


L (m) Cylinder length (not including end caps) 


(m) 


Cylinder cavity length 


Lp (m) 


Floor width (wall to wall) 


Circumferential grid coordinate (non-dimensional) 
for point (k, H) 


M 


Number of longitudinal half-wavelengths for 
structural mode 


m = pt (kg/m^) Average surface mass of the shell or floor 


iHg (kg/m ) Mass/unit area of end cap trim lining, facing the 
cavity 

mp (kg/m ) Mass/unit area of cylinder wall trim lining, 
facing the cavity 


N 


Propeller rpm 


N 


Structural mode number counter only, associated 
with mode r = (M,N) 


LIST OP SYMBOLS 
(Continued) 

NI Number of one-third octave band center frequencies 

(Max 26) 

n Acoustic mode number n = (q,i) 

n Number of shell circumferential waves 

s 

n^ Number of half -wavelengths In floor width 

n* Number of terms in displacement series for shell 

or floor 

p[(k,il),t] (Pa) Propeller pressure at time t at grid location 

(k, £) 

q Number of longitudinal half-wavelengths for 

acoustic mode, associated with mode n = (q,l) 

Rp (m) Propeller radius 

r Structural mode number r = (M,N) 

Tp (m) Radial distance from center of fuselage cylinder 

to the axis of rotation of the propeller 

= BPP ^(secs) Period of propeller noise signature 

t (secs) Time coordinate, over propeller noise period t=0 
to Ti 

t (m) Equivalent thickness of shell or floor (see page 

10 ) 

t^(m) Actual thickness of shell or floor 
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LIST OP SYMBOLS 
(Continued) 


u,v,w Shell and floor degrees of freedom 

V (m^) Cavity volume, above the floor 

I W I (mks rayls) Amplitude of wave impedance of trim insulation 
x,y,z Co-ordinate system (see Figure 3) 


X 


X X 

123 


Local coordinate system used for ANOPP propeller 
noise prediction (see Figure 3) 


Zp (m) Location of propeller relative to the front of the 

cylinder (of length L) 


a(dB/m) Trim insulation attenuation constant 

A = Grid spacing for propeller noise predictions 

Acoustic mode n normalization parameter (see 

Ref. [2]. 

End cap trim lining loss factor (in flexure) 

Acoustic loss factor for mode n 

Average one-third octave band acoustic loss 
factor 


Cylinder trim lining loss factor (in flexure) 

^struc Average one-third octave band structural loss 
factor 


^rad 

r 


Average one-third octave band external radiation 
loss factor 
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LIST OP SYMBOLS 
(Continued) 

Total structural loss factor. Including effect of 
trim 

Average one-third octave band total structural 
loss factor, including trim 

hp’ Internal radiation loss factor, due to closely 

coupled structural and acoustic modes (see Page 
^3, Ref. [2]). 

0( degrees) Shell angular coordinate (see Figure 3)* 6 = 0 is 

at fuselage bottom centerline 

0Q(degrees) Angle at which shell/floor joint is located 

Acoustic wavelength in trim insulation 

Poisson's ratio for shell or floor 

V Poisson's ratio for end caps 

p(lcg/m^ ) Density of shell or floor material 


Trim transmission coefficient = -10 log (t,) dB 

4> (degrees) Angular position of propeller hub relative to 
fuselage bottom centerline (see Figure 2) 

<P (degrees) Phase of wave impedance of trim insulation 

. ( k , «. ) 

? pj (degrees) Phase of Fourier component of propeller 

pressure signature at propeller harmonic H and 
grid location (k,il) 


> 1 ^ 3 ( 2 , 0) Shell displacement at location (z, 0) in mode r 
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LIST OP SYMBOLS 
(Continued) 


\|^p(z,x) Floor displacement at location (z,x) In mode r 

4'„(r,H) Generalized modal forcing function due to 

propeller noise, mode r at propeller harmonic H 

0) (radians/sec) Angular frequency 

( radlans/sec ) Angular frequency of propeller harmonic H 
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